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ABSTRACT 

Understanding  the  heating  of  electrons  to  quasi-thennal  energies  at  collisionless  shocks  has  broad  implications  for 
plasma  astrophysics.  It  directly  impacts  the  interpretation  of  X-ray  spectra  from  shocks,  is  important  for  understanding 
how  energy  is  partitioned  between  the  thermal  and  cosmic-ray  populations,  and  provides  insight  into  the  structure  of  the 
shock  itself.  In  previous  work  by  Ghavamian  et  al.  we  presented  observational  evidence  for  an  inverse-square  relation 
between  the  electron-to-proton  temperature  ratio  and  the  shock  speed  at  the  outer  blast  waves  of  supernova  remnants  in 
partially  neutral  interstellar  gas.  There  we  outlined  how  lower  hybrid  waves  generated  in  the  cosmic-ray  precursor  could 
produce  such  a  relationship  by  heating  the  electrons  to  a  common  temperature  independent  of  both  shock  speed  and  the 
strength  of  the  ambient  magnetic  field.  Here  we  explore  the  mechanism  of  lower  hybrid  wave  heating  of  electrons  in 
more  detail.  Specifically,  we  examine  the  growth  rate  of  the  lower  hybrid  waves  for  both  the  kinetic  (resonant)  and 
reactive  cases.  We  find  that  only  the  kinetic  case  exhibits  a  growing  mode.  At  low  Alfven  Mach  numbers  the  growth  of 
lower  hybrid  waves  can  be  faster  than  the  magnetic  field  amplification  by  modified  Alfven  waves. 

Subject  headings:  cosmic  rays  —  shock  waves  —  supernova  remnants 


1.  INTRODUCTION 

The  main  accelerators  of  cosmic  rays  (CRs)  are  widely  be¬ 
lieved  to  be  high-Mach  number  shocks  in  collisionless  plasma, 
here  loosely  defined  as  plasma  where  charged  particles  interact 
predominantly  through  plasma  waves  rather  than  by  Coulomb 
collisions  (for  a  thorough  conceptual  and  historical  review,  see 
Maikov  &  Drury  2001).  However,  a  consensus  is  emerging  that 
CRs  are  not  simply  a  by-product  of  collisionless  shocks,  but  in 
fact  play  an  integral  role  in  the  shock  structure,  dynamics,  and 
energetics.  For  example,  sound  waves  in  a  CR  pressure  gradient 
can  smooth  out  the  shockjump  in  CR-modified  shocks  (Drury  & 
Falle  1 986).  More  recent  analytic  and  numerical  work  has  shown 
that  modified  Alfven  waves  in  the  CR  precursor  may  amplify  the 
magnetic  field  to  many  times  its  ambient  value  by  generating 
perpendicular  magnetic  field  from  an  initially  quasi-parallel  ge¬ 
ometry  (Lucek  &  Bell  2000;  Bell  &  Lucek  2001;  Bell  2004, 
2005).  Observational  support  for  dramatic  magnetic  field  am¬ 
plification  ahead  of  shocks  exists  in  the  form  of  extremely  thin 
X-ray  synchrotron  rims  of  supernova  remnants  (SNRs)  such 
as  Cassiopeia  A  (Vink  &  Laming  2003),  SN  1006  (Long  et  al. 
2003;  Yamazaki  et  al.  2004),  and  Tycho’s  SNR  (Warren  et  al. 
2005;  Cassam-Chena'iet  al.  2007).  As  noted  by  Cassam-Chenal 
et  al.  (2007),  there  are  two  possible  interpretations  for  the  narrow 
width;  however,  both  require  a  dramatically  amplified  magnetic 
field  ahead  of  the  shock.  Either  the  rims  are  thin  because  the  high 
magnetic  field  causes  rapid  synchrotron  cooling  of  the  X-ray- 
emitting  electrons  or  the  scale  of  the  rims  represents  the  scale  of 
magnetic  field  de-amplification  behind  the  shock. 
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In  this  paper  we  explore  another  area  where  CRs  may  influence 
the  properties  of  the  shock,  namely,  through  the  heating  of  quasi- 
thermal  electrons.  For  shocks  in  collisionless  plasma  the  heating 
of  electrons  must  occur  through  the  damping  of  waves  generated 
by  the  other  more  massive  charged  particles  that  dominate  the 
energetics.  Given  the  wide  array  of  possible  plasma  instabilities 
at  collisionless  shocks,  an  observational  relationship  for  electron 
temperature,  Te ,  at  the  shock  front  was  required  to  limit  theoretical 
discussions.  An  inverse  relationship  between  the  initial  ratio  of 
electron  to  proton  temperatures  immediately  behind  the  shock, 
( Te/Tp)0 ,  and  the  shock  velocity,  vs ,  has  been  reported  in  a  series 
of  observational  papers  on  SNR  shocks  (Ghavamian  et  al.  2001, 
2002,  2003,  2007;  Rakowski  et  al.  2003;  Rakowski  2005)  and 
has  also  been  noted  among  the  higher  Alfven  Mach  number 
events  in  a  sample  of  solar  wind  shocks  (Schwartz  et  al.  1988). 
In  Ghavamian  et  al.  (2007),  we  focused  on  shocks  propagating 
into  partially  neutral  gas.  Here  the  collisional  excitation  of  broad 
and  narrow  Balmer  line  emission  at  the  shock  front  can  be  used 
to  diagnose  (Te/Tp)0.  In  Ghavamian  et  al.  (2001,  2002,  2007) 
we  described  the  method  of  simultaneously  constraining  vs  and 
(Te/Tp)0  via  measurement  of  the  width  of  the  broad  Balmer  line 
and  the  ratio  of  broad  to  narrow  Balmer  line  flux  (see  also  Heng 
&  McCray  2007;  Heng  et  al.  2007).  Our  results  are  consistent 
with  an  inverse-square  relationship,  (Te/Tp)0  oc  \!v2s,  for  shock 
speeds  above  ~400  km  s-1  (Ghavamian  et  al.  2007).  Given  that 
Tp  oc  t>2  at  the  shock  front  by  the  Rankine-Hugoniot  jump  con¬ 
ditions,  the  inverse  relationship  between  equilibration  and  shock 
speed  implies  that  the  electron  temperature  itself  is  nearly  constant, 
~0.3  keV,  independent  of  shock  speed. 

The  insensitivity  of  electron  temperature  to  shock  velocity 
suggests  a  heating  mechanism  within  the  extended  diffusive  CR 
precursor  ahead  of  the  shock.  In  this  case  the  electron  heating 
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would  be  more  reflective  of  the  generic  properties  of  CR  accel¬ 
eration  and  diffusion  than  tied  to  the  specific  attributes  of  the 
shock.  In  contrast,  prior  work  on  heating  by  shock-reflected  ions 
that  are  confined  to  within  a  gyroradius  of  the  (quasi-perpendicular) 
shock  (Cargill  &  Papadopoulos  1988)  suggested  that  ( Te/Tp)0 
would  remain  constant  with  shock  velocity.  In  Ghavamian  et  al. 
(2007)  we  suggested  that  lower  hybrid  waves  in  the  CR  precursor 
of  a  perpendicular  shock  might  be  a  plausible  electron  heating 
mechanism. 

Lower  hybrid  waves  are  electrostatic  ion  waves  directed  nearly 
perpendicular  to  the  magnetic  field  with  a  frequency  equal  to  the 
geometric  mean  of  the  electron  and  ion  gyrofrequencies.  Elec¬ 
trons  that  would  otherwise  screen  the  ion  oscillation  are  pinned 
to  magnetic  field  lines.  In  addition,  the  group  velocity  parallel  to 
the  magnetic  field  greatly  exceeds  the  group  velocity  perpendicular 
to  the  field  (u>/k y  :»  u>/k±).  Therefore,  the  wave  can  simulta¬ 
neously  resonate  with  ions  moving  across  the  field  lines  and 
electrons  moving  along  the  field  lines,  facilitating  collisionless 
energy  exchange  between  them.  Based  on  simple  arguments  about 
the  width  of  the  CR  precursor  and  the  electron  diffusion  along  the 
field  lines,  we  showed  that  electron  heating  from  lower  hybrid 
waves  in  the  CR  precursor  would  be  independent  of  both  the  shock 
speed  and  the  magnetic  field.  Here,  we  explore  this  mechanism 
in  more  physical  detail. 

In  §  2  we  calculate  the  growth  rate  of  lower  hybrid  waves,  first 
examining  the  kinetic  (resonant)  case  then  the  reactive  (non¬ 
resonant)  case.  The  treatment  here  is  mathematically  similar  to 
the  work  on  modified  Alfven  waves  by  Achterberg  (1983)  and 
Bell  (2004,  2005)  involving  the  CR  contribution  to  the  plasma 
dielectric  tensor.  The  analysis  also  draws  on  the  work  of  Laming 
(2001a,  2001b)  on  lower  hybrid  waves  from  shock-reflected  ions. 
We  compare  these  growth  rates  with  those  for  magnetic  field  am¬ 
plification  to  assess  the  conditions  under  which  electron  heating 
might  occur.  In  §  3  we  discuss  the  structure  of  the  CR  shock  pre¬ 
cursor  in  more  detail.  We  pay  particular  attention  to  the  magnetic 
field  geometry,  since  the  excitation  of  lower  hybrid  waves  requires 
a  quasi-perpendicular  shock.  We  show  schematically  how  mag¬ 
netic  field  amplification  and  lower  hybrid  wave  heating  might 
coexist  in  the  shock  precursor  for  either  parallel  or  perpendicular 
initial  geometries.  We  also  review  some  other  ideas  for  electron 
heating  and  make  some  quantitative  predictions  from  our  model 
for  various  shock  parameters.  Included  in  the  appendices  are  a 
discussion  of  CR  diffusion  coefficients  and  a  derivation  of  the 
resonant  growth  rate  for  electromagnetic  waves. 

2.  COSMIC-RAY  GROWTH  RATE 
OF  LOWER  HYBRID  WAVES 

Lower  hybrid  waves  ahead  of  collisionless  shocks  have  par¬ 
ticularly  interesting  properties.  They  can  have  a  group  velocity  away 
from  the  shock  equal  to  the  shock  velocity  itself  ( McClements  et  al. 
1997).  This  can  in  principle  allow  the  waves  to  grow  to  large  am¬ 
plitudes,  even  if  their  intrinsic  growth  rate  is  small.  To  determine 
if  lower  hybrid  waves  can  heat  the  electrons  to  the  ~0.3  keV 
temperature  observed,  we  must  first  calculate  the  growth  rate  of 
this  instability  to  see  if  it  will  have  sufficient  power  to  overcome 
the  damping  effect  of  the  electrons  as  well  as  to  compete  with 
other  instabilities  in  the  precursor.  We  calculate  this  growth  rate 
in  both  kinetic  and  reactive  limits,  i.e.,  either  considering  the  CRs 
with  energies  in  resonance  with  the  lower  hybrid  wave  frequency 
or  the  integrated  contribution  of  the  entire  distribution,  respec¬ 
tively  (see,  e.g.,  Melrose  1986).  Related  kinetic  and  reactive  cases 
were  calculated  in  Laming  (2001a,  2001b)  but  only  for  the  case  of 
shock-reflected,  nonrelativistic  ions  gyrating  around  the  magnetic 


field,  represented  as  a  particle  beam.  Here  we  begin  the  discussion 
with  the  resonant  case. 

2.1.  Kinetic  Growth  Rate 

We  model  the  nonnalized  CR  distribution  function  diffusing 
upstream  as 


f(p)  = 


a?/r  (2k  —  3)T(k) 
4\/2(7r/t)3/2/7^  T(k  —  1/2) 


1  + 


(px-mvsy  +pl+pl 


2  Kp2 


(1) 


where  the  coordinate  system  is  aligned  so  that  the  shock  speed  vs 
lies  in  the  x-direction.  Here,  p,  is  defined  as  the  “thermal”  mo¬ 
mentum,  which  we  take  to  be  (3/4 )mvs,  and  n'CR  denotes  the 
density  of  suprathermal  particles  with  distribution  function  f(p). 
The  functional  form  above,  known  as  a  “kappa”  distribution, 
is  often  seen  for  particle  distribution  functions  associated  with 
shocks  in  the  solar  wind  and  may  be  derived  as  equilibrium  dis¬ 
tributions  for  a  system  of  particles  and  waves  under  certain 
conditions  (see,  e.g.,  Laming  &  Lepri  2007),  in  contrast  to  a  sys¬ 
tem  of  particles  which  only  gives  a  Maxwellian  distribution  of 
width  p,.  The  kappa  distribution  resembles  a  Maxwellian  for 
P  <Pt  and,  in  fact,  for  k  — >  oo  is  exactly  a  Maxwellian.  At 
higher  particle  momenta  it  tends  smoothly  to  a  distribution/  oc 
P~2k.  Below,  we  shall  take  At  =  2  to  model  the  well-known 
f(p)  oc  p~4  CR  distribution  predicted  by  diffusive  shock  accel¬ 
eration  in  shocks  with  a  compression  ratio  of  4.  In  connecting  the 
CRs  to  the  lower  energy  particles  in  this  way,  we  are  somewhat 
blurring  the  distinction  between  “CRs”  and  other  suprathermal 
particles  reflected  from  the  shock.  Hence,  we  denote  the  com¬ 
bined  density  of  these  particles  as  n'CR  to  distinguish  it  from 
density  of  true  CRs,  ajcr,  which  will  appear  in  expressions  de¬ 
rived  by  other  authors.  Note  that  all  these  particles  are  distinct 
from  the  ambient  thermal  plasma  upstream  of  the  shock,  which 
here  is  considered  to  be  a  Maxwellian  with  much  lower  tem¬ 
perature  than  p,  used  in  equation  (1)  for  the  upstream  supra- 
thermals.  When  discussing  the  kinetic  instability  we  focus  on 
particles  that  obey  a  diffusion  equation  ahead  of  the  shock,  due 
to  their  interaction  with  turbulence,  rather  than  gyrate  around 
field  lines.  We  qualify  this  distinction  more  carefully  below  in 
our  discussion  of  the  reactive  instability. 

The  appropriate  dispersion  relation  for  equation  (1)  can  be 
found  from  the  cold  plasma  dielectric  tensor.  For  electrostatic 
waves  at  frequencies  close  to  the  lower  hybrid  wave  frequency 
we  have  (Laming  2001a,  2001b) 


->  2  2 

LO  „  UJ  •  U) 

Pe  •  2/q  pi  pe n2c 


Kl  =  1  +  -/f  sin26>  -  -f-  -  cos 


K 

Anq1 


id- 

k  •  df  /dp 
u>  —  k  •  v 


dJp  =  0, 


(2) 


where  9  denotes  the  angle  between  the  wavenumber  of  the  per¬ 
turbation  and  the  preshock  magnetic  field \u)px  =  (4nq2nx/mxy1/2> 
is  the  plasma  frequency  of  a  given  species  x  (electrons,  ions, 
CRs,  etc.)  with  charge  q,  density  nx,  and  mass  mx;  Qx  =  qBI 
(aw7Lc)  is  the  cyclotron  frequency  of  species  x  (with  7l  being 
the  Lorentz  factor);  and  unadorned  w  is  the  lower  hybrid  wave 
frequency  which  is  the  geometric  mean  of  the  electron  and  pro¬ 
ton  cyclotron  frequencies.  Using  the  Landau  prescription  for 
evaluating  the  integral  at  the  resonant  pole  and  taking  only  the 
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imaginary  parts  of  the  dielectric  tensor  equation,  we  find  the 
growth  rate  for  the  lower  hybrid  waves, 


7  = 


2ttV  w2 
k2  u)2pi  +  u2pe  cos2# 


k  •  v)k  • 


(3) 


1986;  Pelletier  et  al.  2006) 

37T  kli  77CR  /  COS  #  4  7T  \ 

7-S, res  =  77- - 1 - tttLs  cos  v  -  -  -jvs  sm(p)k\\rg, 

16  va  n i  \|  cos  #|  3  4  /  " 

(8) 


To  compute  7,  we  take  n'CR  oc  exp (—xvs/D),  where  D  is  the 
CR  diffusion  coefficient  and  is  in  principle  dependent  on  the  CR 
momentum  (making  /  =  D/vs  the  characteristic  diffusive  length 
scale).  With  this  substitution  we  begin  the  evaluation  of  the  last 
integral  in  equation  (3), 


/” 


iTTf^fdp^  = 


lCR 


(2k  —  3)T(k)  27 TK 


4\/2(7rK)3/2/?2  T(k  —  1/2)  k  —  1 


xp~t 


1  + 


(Px  ~  mvs) 


2“i  1  ~K 


2kP2 


0-xvs/D 


(4) 


where  v\  is  the  Alfven  speed,  n is  the  density  of  ions  in  the  back¬ 
ground  plasma,  and  rg  is  the  gyroradius  of  CRs.  This  expression 
differs  from  that  in  the  cited  references  in  the  factor  cos  and  the 
term  in  sin  tj>,  where  <f>  is  the  angle  between  the  shock  velocity 
and  the  magnetic  field.  At  perpendicular  shocks,  the  growth  rate 
of  resonant  Alfven  waves  can  be  neglected,  but  at  parallel  shocks, 
may  be  larger  than  that  for  lower  hybrid  waves,  depending  on  the 
ratio  «£R/«cr.  However  as  we  shall  argue  below,  all  shocks  sub¬ 
ject  to  magnetic  field  amplification  become  perpendicular,  and 
this  is  the  geometry  where  lower  hybrid  waves  are  most  effectively 
excited,  so  we  neglect  from  here  onward.  Bell  (2004)  dis¬ 
covered  a  nonresonant  growth  rate  for  Alfven  waves,  with  ap¬ 
proximate  growth  rate 


where  we  have  separated  out  the  components  of  p  perpendicular 
to  the  shock  (and  k)  from  px.  Substituting  back  into  equation  (3) 
yields  a  growth  rate 
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For  waves  to  stay  in  contact  with  the  shock,  uj/k  ~  —2vs(  Laming 
2001a)  in  the  cold  plasma  electrostatic  limit,  i.e.,  reflected  ions 
returning  to  the  shock  excite  the  waves.  This  remains  generally 
true  when  these  approximations  are  relaxed  (Laming  2001b), 
so  the  /-function  picks  out  px  =  —2 mvs  (CRs  returning  to  the 
shock).  Hence, 


Ti¬ 

ll 

'7T\3/2 

^  P,( 
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Assuming  dD/dpx  =  0  (see  Appendix  A  on  CR  diffusion  co¬ 
efficients)  and  k  =  2,  only  the  first  tenn  within  the  square  brack¬ 
ets  remains.  Evaluating  the  T  functions  and  substituting  pt  = 
(3/4)77777,  we  arrive  at  the  following  expression  for  the  kinetic 
growth  rate  of  lower  hybrid  waves, 
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where  uj'pi  denotes  the  plasma  frequency  for  n'CR.  Substituting  in 
the  frequency  definitions  we  note  that  approximately  7  cx  (t7^r/ 
m)kli.  Before  proceeding,  we  pause  to  compare  this  growth  rate 
with  those  for  magnetic  field  amplification.  In  the  case  where 
Alfven  waves  are  resonantly  excited,  the  growth  rate  is  (Melrose 


which  has  a  maximum  value  of  MaIX/zcr^t/,.  According  to 
Bell  (2005),  this  instability  operates  for  arbitrary  orientations  of 
B ,  vs,  and  fc,  indicating  that  it  will  also  amplify  magnetic  field  at 
perpendicular  shocks.  Its  growth  rate  is  strongest  for  fc|  \B  and 
zero  for  k  _L  B.  Equating  the  maximum  value  of  7_a:„onres  with 
the  lower  hybrid  wave  growth  rate  calculated  above,  we  find  the 
critical  Alfvenic  Mach  number  Ma  —  l2n'CR/ncR,  such  that  for 
higher  Ma,  CRs  preferentially  amplify  magnetic  field,  and  for 
lower  Ma,  they  generate  lower  hybrid  waves.  The  numerical  value 
depends  on  the  ratio  n'CR/ncR.  In  §  2.2  we  argue  that  these  two 
densities  should  not  be  the  same  and  that  77 /.R  >  77  cr,  following 
from  a  consideration  of  the  reactive  growth  rate  for  lower  hybrid 
waves. 


2.2.  Reactive  Growth  Rate 

The  reactive  case  involves  the  integrated  contribution  to  the 
growth  rate  from  the  entire  CR  distribution.  Thus,  we  examine 
successive  orders  in  an  expansion  of  f(p)  to  see  if  they  produce 
any  growing  modes.  Although  no  instability  is  found  in  this  pro¬ 
cess,  we  do  uncover  potentially  interesting  constraints  on  the 
properties  of  f(p). 

We  consider  again  the  last  term  in  equation  (2),  the  CR  con¬ 
tribution  to  the  dielectric  tensor  (e.g.,  Melrose  1986),  which  in¬ 
cludes  a  factor  that  reduced  to  unity  for  the  resonant  case, 
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For  the  case  of  a  beam  of  CR  particles  localized  around  vs,  one 
recovers  the  usual  beam  reactive  instability  (Laming  2001b, 
eq.  [A6]).  However,  as  we  demonstrate  below,  for  a  more  phys¬ 
ical  quasi-isotropic  CR  distribution  no  instability  is  recovered. 
We  consider  the  case  where  CRs  drift  with  velocity  vs  and  waves 
are  generated  with  k  ||  vs.  The  CR  distribution  function  from 
Appendix  A,  expanded  in  tenns  of  the  cosine  of  the  angle  be¬ 
tween  k  and  u  or  k  and  vs,  cos  a  =  p,,  is 
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such  that  f  pvf  dkl  =  fovs  as  before.  Neglecting  terms  of  order 
v\ ]/v2,  the  CR  contribution  to  the  dielectric  tensor  becomes 
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Expanding  the  numerator  into  terms  divisible  by  (oj/kv  —  p)  and 
evaluating  the  integral  over  p,  we  obtain 
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We  then  evaluate  this  integral  in  the  two  limiting  cases  away 
from  the  pole,  lo^>  kv  and  ui  <SC  kv.  For  u  kv. 
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All  terns  in  brackets  cancel  save  for  one,  giving 
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The  first  term  goes  to  zero  so  long  as /o(oo)  — »  0  faster  than p  3, 
and  the  second  term  is  otncRlfi/m),  leaving 
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This  simply  modifies  the  1  luf  tenn  in  the  dielectric  tensor,  chang¬ 
ing  the  frequency  of  the  solution  but  not  creating  any  complex 
roots;  hence,  no  instability  is  generated. 

Likewise,  in  the  case  where  lo  <SC  kv. 
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All  the  terms  but  the  first  are  negligible  in  this  limit;  hence, 
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This  leads  to  a  full  dispersion  relation  that  can  be  written  as 
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which  also  lacks  complex  roots,  regardless  of  the  actual  evalu¬ 
ation  of  the  integral. 

Higher  order  tenns  in  the  expansion  of f  =  fo(  1  +  p(3vs/v)  + 
p2(3vs/v)2 12  +  p3(3vs/v)2 16  +  ■  ■  .)  give  rise  to  higher  order 
terms  in  lo.  For  lo  kv,  the  dispersion  relation  equation  (18)  be¬ 
comes  to  next  highest  order 


(23) 


which  is  stable  since  the  terms  oc J(l/v)(dfcR/dp)p2  dp  are  of 
order  ~w2c Rv2/v2  <SC  lo2c r.  When  lo  <3C  kv,  the  dispersion  rela¬ 
tion  equation  (22)  takes  on  the  next  highest  order  terms 
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to  become  a  quartic  equation.  This  has  four  real  solutions  so  long 
as  f(l/v5)(dfo/dp)p2  dp  >  0. In  fact  if  J(l/v5)(dfo/dp)p2  dp  < 
0,  the  addition  of  higher  order  tenns  in  the  expansion  of  the  CR 
distribution  function  would  dramatically  alter  the  character  of 
the  solutions,  a  situation  that  must  be  considered  unphysical.  We 
require  J(l/v5)(dfo/dp)p2  dp  >  0,  which  means  that  at  low  mo¬ 
menta,  dfo  I  dp  >  0,  and  the  CR  distribution  cannot  be  monoton- 
ically  decreasing  from  v  =  p  =  0.  Our  use  of  the  kappa  distribution 
in  §  2. 1  may  therefore  be  questioned.  However  the  resonance  at 
px  =  —2mvs  places  it  well  into  the  region  of  the  distribution 
where  the  gradient  is  negative,  and  so  modifications  to  the  low- 
momentum  behavior  would  have  very  little  effect  on  our  result. 
However,  this  observation  does  imply  that  the  distribution  of 
particles  obeying  a  diflusion  equation  ahead  of  the  shock  is  un¬ 
likely  to  extend  down  to  zero  momentum.  Some  natural  break 
must  exist  between  the  quasi-thermal  population  gyrating  around 
field  lines  and  the  CRs  diffusing  in  turbulence.  The  foregoing 
work  also  neglects  the  CR-induced  current  in  the  background 
plasma.  The  inclusion  of  such  effects  leads  to  the  modification 
w2;-  — >  ui2pi  —  uipCR  and  has  no  effect  on  reactive  instabilities. 

3.  DISCUSSION 

3.1.  Electron  Heating  or  Magnetic  Field  Amplification? 

We  have  calculated  the  growth  rate  for  waves  that  damp  by 
heating  electrons,  in  a  CR  shock  precursor  using  similar  ap¬ 
proximations  and  techniques  to  those  employed  by  Bell  (2004). 
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Fig.  1 . — Schematic  illustrations  of  the  amplification  of  magnetic  field  by  the 
nonresonant  modified  Alfven  waves  in  the  shock  precursor  in  the  parallel  orientation 
of  the  ambient  field  with  respect  to  the  shock  normal.  The  evolution  of  a  single  field 
line  in  an  exponential  purely  growing  mode  is  shown.  As  the  field  is  amplified,  the 
shock  becomes  quasi-perpendicular  and  the  effective  Ma  decreases,  eventually  to 
the  point  where  the  lower  hybrid  wave  growth  takes  over,  allowing  a  short  region 
of  electron  heating. 


Both  the  lower  hybrid  wave  heating  of  electrons  and  the  growth 
of  magnetic  field  through  modified  Alfven  waves  redistribute 
energy  within  the  CR  precursor.  An  important  question  is  which 
of  these  is  more  effective,  i.e.,  which  grows  faster?  Above,  we 
derived  a  critical  Alfven  Mach  number,  Ma  —  \2n'CR/ncR,  which 
divides  the  regime  of  magnetic  field  amplification  from  that  of 
lower  hybrid  wave  growth.  Following  from  the^treatment  of  the 
reactive  instability  above,  we  estimate  «cr  =  Sp^fcvA'Kp2  dp  ~ 
{6hT)n'CRvs/vinj,  where /cr  is  given  by  equation  ( 1 )  and p ;nj  is  the 
injection  momentum  where  particles  may  begin  to  participate  in 
a  diffusive  shock  acceleration  process.  The  approximate  result 
nc r/”cr  —  Viaj/2vs  gives  MA  ~  &vm^/vs  as  the  critical  Alfven 
Mach  number. 

The  next  step  in  determining  the  critical  Alfven  Mach  number 
is  to  find  an  appropriate  ujnj  for  the  injection  of  seed  particles  into 
the  CR  acceleration  process.  Zank  et  al.  (2006)  argue  that  quasi¬ 
perpendicular  shocks  have  similar  injection  requirements  to 
quasi-parallel  shocks  (ujnj  ~  2vs),  but  that  highly  perpendicular 
shocks  require  much  higher  injection  energies.  In  the  case  of 
nonresonant  magnetic  field  generation,  we  also  consider  the  case 
of  a  highly  perpendicular  shock,  since  the  generated  magnetic 
field  will  be  perpendicular  and  much  stronger  than  the  initial  mag¬ 
netic  field.  Zank  et  al.  (2006)  give  the  injection  velocity  as 
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where  r  is  the  shock  compression  ratio,  and  Z>Bohm  and  D±  are 
the  CR  diffusion  coefficients  in  the  Bohm  limit  and  in  the  perpen¬ 
dicular  direction,  respectively.  Reville  et  al.  (2008)  giveZ)Bohm/ 
D±  ~  3  for  CRs  where  krg  ~  1,  and  so  Djnj  ~  10ds.  Thus,  the 
Alfven  Mach  number  at  which  lower  hybrid  wave  growth  takes 
over  from  magnetic  field  amplification  should  be  about  12-60, 
unless  the  magnetic  field  saturates  at  a  lower  value  (i.e.,  higher 
Ma)  before  this  is  reached. 

The  growth  of  lower  hybrid  waves  is  most  efficient  at  a  quasi¬ 
perpendicular  shock,  whereas  the  growth  of  magnetic  field  through 


Fig.  2. — Same  as  Fig.  1,  but  for  an  initially  perpendicular  shock.  The  evolution 
of  a  purely  growing  mode  is  illustrated.  The  magnetic  field  amplification  is  less 
strong  than  in  the  quasi-parallel  case,  and  the  shock  geometry  remains  quasi¬ 
perpendicular. 


modified  Alfven  waves  is  strongest  at  a  quasi-parallel  shock  (Bell 
2005).  This  apparent  contradiction  is  actually  easily  resolved.  At 
an  initially  quasi-parallel  shock,  Bell  (2005),  Reville  et  al.  (2008), 
and  Zirakashvili  et  al.  (2008)  show  that  a  highly  helical  mag¬ 
netic  field  develops.  The  distortion  of  an  initially  parallel  field 
line  is  shown  schematically  in  Figure  1,  showing  the  evolution  of 
the  shock  from  quasi-parallel  to  quasi-perpendicular.  A  similar 
schematic  in  Figure  2  shows  the  evolution  of  an  initially  quasi¬ 
perpendicular  shock,  where  magnetic  field  is  amplified  orthogonal 
to  the  preexisting  magnetic  field,  but  where  the  shock  remains 
quasi-perpendicular.  In  both  cases  a  perpendicular  field  is  gener¬ 
ated,  thus  allowing  lower  hybrid  wave  growth  and  electron  heat¬ 
ing  in  a  region  close  to  the  shock  as  indicated. 

Another  potential  problem  is  the  cavities  seen  in  simulations 
of  the  growth  of  modified  Alfven  waves  (e.g.,  Bell  2005).  The 
helical  field  from  an  initially  quasi-parallel  geometry  naturally 
creates  a  filamentary  structure,  dragging  the  thermal  plasma  with 
it,  while  CRs  tend  to  accumulate  in  the  low-density  cavities.  This 
is  problematic  for  our  mechanism  that  requires  spatial  coinci¬ 
dence  between  CRs,  magnetic  field,  and  thermal  plasma.  A 
possible  solution  is  that  the  growth  of  lower  hybrid  waves  takes 
over  from  the  growth  of  modified  Alfven  waves,  so  that  the  CR- 
driven  magnetic  field  never  reaches  its  final  saturated  state.  Bell 
(2004),  Reville  et  al.  (2008),  and  Zirakashvili  et  al.  (2008)  derive 
a  saturation  magnetic  field  by  setting  7s,nonres  =  0  in  equation  (9) 
to  give  SB  ~  jrg/4iv  or  SB2/ 87t  ~  «cR"TAAinj/2.  This  gives  an 
Alfven  Mach  number  at  saturation  of  M2k  ~  n,/10ncR  (assuming 
Djnj  ~  10O,  which  for  likely  parameters  »,/« cr  ~  103  gives  a 
value  of  Ma  of  similar  magnitude  to  but  possibly  lower  than  that 
where  the  electron  heating  is  expected  to  take  over.  Bearing  in 
mind  that  we  took  the  strongest  growth  rate  for  magnetic  field 
amplification  to  estimate  where  electron  heating  takes  over,  it  is 
quite  plausible  that  the  amplified  magnetic  field  never  reaches 
saturation.  Also,  as  the  initial  shock  state  becomes  more  quasi¬ 
perpendicular,  the  growth  rate  slows  down,  and  the  circularly 
polarized  Alfven  waves  become  elliptically  polarized,  ultimately 
becoming  linearly  polarized  in  the  limit  of  a  true  perpendicular 
shock,  eliminating  the  growth  of  such  cavities. 

Pelletier  et  al.  (2006)  find  that  the  nonresonant  instability  of 
Bell  (2004,  2005)  dominates  over  the  more  familiar  resonant 
instability  when  the  shock  velocity  vs  is  greater  than  a  few  times 
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ecrC,  where  £cr  is  the  ratio  of  the  CR  energy  density  to  the  ki¬ 
netic  energy  density  of  the  shock.  Niemiec  et  al.  (2008)  simulated 
the  CR-driven  amplification  of  magnetic  field  in  a  parallel  shock 
using  particle-in-cell  simulations,  which  can  naturally  account  for 
the  back-reaction  of  the  generated  magnetic  field  on  the  CR  cur¬ 
rent.  In  conditions  where  the  nonresonant  mode  should  grow,  they 
find  magnetic  field  amplification  only  to  SB  ~  B.  The  magnetic 
field  again  produces  filaments,  but  they  do  not  find  CR  accumu¬ 
lation  in  the  filament  cavities.  They  do  not  find  strong  growth  and 
argue  that  saturation  occurs  because  the  incoming  flow  to  the 
shock  is  decelerated  by  the  CRs,  reducing  their  relative  velocity 
and,  hence,  the  CR  current. 

For  our  electron  heating  model,  the  precise  degree  of  magnetic 
field  amplification  is  unimportant  so  long  as  the  CR  diffusion 
coefficient  remains  proportional  to  1  IB.  It  is  only  necessary  that 
the  shock  be  sufficiently  quasi-perpendicular  to  allow  CRs  to 
generate  lower  hybrid  waves.  A  reduced  CR  current  does  not 
necessarily  produce  an  appreciable  affect  on  the  kinetic  growth 
rate  for  lower  hybrid  waves.  So  long  as  the  current  does  not  van¬ 
ish,  the  initial  effect  of  reducing  vs  in  equation  (6)  is  to  bring  more 
CRs  into  resonance  with  the  lower  hybrid  waves.  Another  esti¬ 
mate  of  the  CR  density  necessary  to  heat  electrons  may  come  from 
the  long-wavelength  limit  of  the  magnetic  field  amplification,  when 
7a  =  ( k  •  BncRqvs/riimi)12,  for  both  parallel  and  perpendicular 
cases  (Bell  2005).  Electron  heating  then  requires  7  ~  0,«cr/«,  > 
(>icrvs  cos  <p/riiv i„j) 1/2  H;,  taking  k  =  fii/umj  (probably  an  over¬ 
estimate),  yielding  ncR^i  >  cos  0/10.  At  cos  <f>  <  (me/mi)12, 
the  values  typical  for  lower  hybrid  wave  propagation,  the  value  for 
mcr/h,  is  low  enough  (0.001-0.01)  to  make  electron  heating  by 
CRs  plausible. 


shock  is  turbulent,  the  importance  of  such  electric  fields  is  less 
clear.  Electron  E  x  B  drift  along  the  shock  front  will  result  in 
periods  of  energy  loss  as  well  as  energy  gain  by  the  cross  shock 
potential  and,  hence,  no  net  heating.  It  has  been  argued  (Gedalin 
et  al.  2008)  that  in  certain  cases  the  shock  front  may  be  suffi¬ 
ciently  thin  (  length  scales  of  order  c/u>pe)  that  the  electrons  are 
effectively  demagnetized.  One  might  expect  to  see  electron 
heating  increase  with  shock  velocity  (or  MA)  once  this  condition 
becomes  satisfied.  Examination  of  solar  wind  shocks  suggests 
that  such  thin  shocks  are  rare  at  best  and  certainly  not  ubiquitous. 

We  see  no  evidence  for  an  increase  in  electron  heating  in 
SNRs  up  to  shock  velocities  of 6000  km  s  “ 1  ( 1 E  0 1 02.2— 7219; 
Hughes  etal.  2000),  and  possibly  up  to  20,000  km  s-1  (SN  1993J; 
Fransson  et  al.  1996).  At  higher  Mach  numbers  such  as  those 
expected  in  gamma-ray  burst  afterglows,  the  convective  electron 
gyroradius  may  easily  reduce  to  less  than  the  electron  inertial 
length,  making  the  cross  shock  potential  a  candidate  electron 
heating  mechanism. 

Schwartz  et  al.  (1988)  have  made  a  survey  of  a  number  of 
solar  wind  shock  crossings  observed  in  situ.  They  find  7V7)  ex 
1/Ma  for  Ma  greater  than  about  2-3.  At  lower  Ma,  there  is  a 
wide  scatter  in  7V7)  about  TeITj  ~  1 .  At  these  slower  shocks,  Te 
correlates  very  well  with  the  change  in  ion  velocity  squared,  sug¬ 
gesting  that  both  are  due  to  the  same  mechanism,  presumably  the 
cross  shock  potential.  The  switch  to  TeITi  cx  1  /Ma  atMA  ~  2-3 
is  possibly  due  to  the  onset  of  turbulent  shock  structure  at  higher 
Mach  numbers. 

We  can  explore  the  conditions  required  for  the  validity  of  the 
laminar  approximation  by  adopting  the  criterion  of  Tidman  & 
Krall  (1971)  for  the  existence  of  a  magnetosonic  soliton. 


3.2.  Other  Electron  Heating  Mechanisms 

Several  other  researchers  have  considered  the  generation  of 
waves  in  a  shock  precursor  as  a  means  of  heating  electrons.  Ohira 
&  Takahara  (2007)  and  Shimada  &  Hoshino  (2000)  have  both 
considered  the  model  of  Cargill  &  Papadopoulos  (1988)  in  more 
detail,  using  particle-in-cell  codes  rather  than  a  hybrid  approach. 
Other  references  (Dieckmann  et  al.  2000;  McClements  et  al. 
2001;  Schmitz  et  al.  2002)  focus  more  on  the  electron  injection 
problem  for  diffusive  shock  acceleration,  rather  than  the  thermal 
electron  temperature,  again  invoking  various  wave  modes  in  a 
reflected-ion  precursor.  Our  principal  departure  from  these  works 
has  been  to  treat  similar  wave  modes  upstream  of  the  shock,  but 
excited  by  CRs  undergoing  diffusive  shock  acceleration  rather 
than  by  quasi-thermal  ions  reflected  from  the  shock.  This  allows 
electron  heating  to  occur  over  a  much  more  extended  upstream 
region  dictated  by  the  CR  diffusion  coefficient,  D ,  rather  than  the 
ion  gyroradius.  In  addition,  expressing  the  thickness  of  this  re¬ 
gion  as  /  ~  D/vs  naturally  results  in  electron  heating  that  is  es¬ 
sentially  independent  of  shock  speed,  as  argued  from  observations 
of  Balmer-dominated  shocks  in  Ghavamian  et  al.  (2007).  On  the 
other  hand,  if  the  thickness  of  the  electron  heating  region  is  com¬ 
parable  to  the  ion  gyroradius,  then  /  tx  vs  and  Te  <x  v2.  This  re¬ 
sults  in  constant  7V7)  with  vs,  contrary  to  what  is  observed. 

A  number  of  other  authors  have  investigated  the  role  of  the 
cross  shock  potential  in  heating  the  electrons.  Inside  the  (quasi¬ 
perpendicular)  shock  ramp,  the  magnetic  field  may  “overshoot,  ” 
i.e.,  increase  to  a  value  greatly  in  excess  of  its  asymptotic  down¬ 
stream  strength  before  decreasing  again.  The  electric  field  arising 
from  the  small  charge  separation  associated  with  this  magnetic 
field  gradient,  E  ~  (9/&)(R2)/(87re«,-),  can  decelerate  ions  and 
accelerate  electrons.  Such  effects  are  known  to  be  important  at 
low-Mach  number  shocks  where  a  laminar  approximation  holds 
(e.g.,  Scudder  et  al.  1986).  At  higher  Mach  numbers,  where  the 


Ml 

Mj-l 


<Mi 


2Mf-4Ms  +  3  +  21nMs 
^  (Mj  —  1  —  21nMJ) 2 


(26) 


The  relationship  between  Ma  and  Ms  predicted  by  this  relation 
is  plotted  in  Figure  3.  The  criterion  above  indicates  that  the  lam¬ 
inar  approximation  breaks  down  at  slightly  lower  Mach  numbers 
in  the  solar  wind  shocks  than  indicated  by  the  behavior  of  TJTt 
in  Schwartz  et  al.  (1988).  Magnetic  field  amplification  by  about 
an  order  of  magnitude  in  SNR  shocks  for  the  400  km  s  ~~ 1  shocks 
observed  in  the  Cygnus  Loop,  where  Ghavamian  et  al.  (2007) 
find  complete  electron-ion  equilibration,  would  bring  MA  down 
to  the  same  range  as  indicated  by  Schwartz  et  al.  (1988),  possibly 
suggesting  that  the  cross  shock  potential  is  at  work  for  the  lower 
velocity  SNR  shocks.  At  the  higher  velocity  shocks  in  the 
Schwartz  et  al.  (1988)  sample,  which  are  all  perpendicular,  an 
empirical  relationship  TeITi  tx  1/MA  emerges.  Such  a  behavior 
can  be  consistent  with  our  model  if  we  make  the  assumption  that 
the  CRs  accompanying  solar  wind  shocks  are  nonrelativistic, 
suprathermal  particles.  Then  the  diffusion  coefficients  take  on  an 
extra  factor  vslc,  assuming  that  the  CR  velocity  is  proportional  to 
the  shock  velocity.  This  extra  power  of  the  shock  velocity  in  the 
diffusion  coefficient  results  in  Te  cx  vs  and  Te!Tj  cx  \!vs. 


3.3.  Heating  versus  Damping  of  Lower  Hybrid  Waves 
and  the  Width  of  the  Precursor 

We  have  discussed  whether  growth  of  lower  hybrid  waves 
may  compete  with  CR-induced  magnetic  field  amplification,  but 
have  not  yet  discussed  whether  this  growth  rate  is  sufficient  to 
balance  the  damping  rate  of  lower  hybrid  waves  by  electrons.  To 
answer  this  question  we  compare  the  electron  heating  rate  from 
diffusive  scattering  off  the  lower  hybrid  waves  with  the  energy 
input  into  the  lower  hybrid  waves  from  the  CR  turbulence.  The 
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by  Coulomb  collisions,  with  characteristic  timescale  10 10 (77 
108  K ya lne  s.  Equating  this  to  1 0 8 /« e  s  yields  a  maximum 
temperature  of  T  ~  108(10_2)2/3~  5  x  106  K.  This  is  very  close 
to  the  temperature  found  in  Ghavamian  et  al.  (2007),  0.3  keV,  or 
3.5  x  106  K.  Put  another  way,  the  temperature  found  in  Ghavamian 
et  al.  (2007)  is  consistent  with  electron  heating  such  that  neutral 
hydrogen  can  survive  to  encounter  the  shock  front  proper.  How¬ 
ever,  CR  precursors  at  the  small  end  of  the  range  considered 
above  (~10n  cm)  would  not  allow  any  significant  electron  col- 
lisional  equilibration  to  occur.  Allowing  for  compressional  heat¬ 
ing  of  the  electrons  as  they  go  through  the  shock,  a  precursor 
electron  temperature  of  order  106  K  requires  a  precursor  length 
of~10  7vs/ne  ~  1015(^/1000  km  s_1)  cm  or  a  minimum  CR  dif¬ 
fusion  coefficient  of  D  ~  1023()7/1000  km  s-1)2  cm2  s-1. 

The  electric  field  in  the  lower  hybrid  waves  will  be  given  by 
the  limit  derived  by  Kamey  (1978), 


Fig.  3. — Allowed  range  of  MA  as  a  function  of  Ms  for  the  existence  of  a 
magnetosonic  soliton,  from  Tidman  &  Krall  (1971).  The  upper  limit  is  given  by 
the  solid  line,  and  the  lower  limit  is  given  by  the  dashed  line. 
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electron  heating  rate  per  unit  area  of  shock  is  nefRmen\\\\tvs/2, 
where fR  ~  exp  (—  w2 12k^v2te y)  ~  exp(— 2v2smilirteme)  is  the  frac¬ 
tion  of  electrons  in  resonance  with  the  lower  hybrid  waves  (using 
uj/k  ~  2vs),  Kim  is  the  parallel  electron  velocity  diffusion  coef¬ 
ficient  in  lower  hybrid  turbulence,  and  t  is  the  period  of  time 
spent  by  an  electron  in  the  turbulence.  This  time  t  =  l/vs,  where 
/  is  the  precursor  depth.  The  electron  heating  is  balanced  by  en¬ 
ergy  input  to  the  turbulence  by  CRs  with  rate  2^/Etm^l.  Putting 
Eturb  =  {8E2 /Sir)2 u!~e/Q2e  and  ac || ||  =  q28E2kp4mlk±u>,  we  de¬ 
duce  a  growth  rate  7  =  fRq2B2 /\6mjineujc2  =  fRuil  16.  The  ki¬ 
netic  growth  rate  derived  earlier  in  equation  (7),  in  units  of  O,-,  is 
proportional  to  «cr/m/.  Thus,  as  long  as  this  ratio  is  comparable 
to  or  larger  than  the  fraction  of  electrons  that  are  in  resonance 
with  the  lower  hybrid  waves  the  growth  rate  outlined  above  will 
be  sufficient  to  heat  the  electrons. 

Another  constraint  on  /zcr/«  ;  comes  from  equation  (9).  For  mag¬ 
netic  field  amplification,  we  require  hcrTz,  >  k\\v\/(vs£li).  Since 
A'||  >  \/rg,  the  gyroradius  of  CRs  at  injection,  «cr/m,-  >  v\/ 
Hnji7  ~  1/10M^  ~  10-3  for  Ma  ~  10  and  vm^  ~  10t>s.  Taking 
the  maximum  growth  rate  estimated  from  equation  (9),  7  = 
MtSlincv.l2rii  ~  O,/20Ma,  where  /,■  =  du)pi  =  va/^i  is  the  ion 
inertial  length,  we  estimate  a  characteristic  length  of  17/7  = 
20M^//~  5  xlO10  cm.  This  requires  a  CR  diffusion  coefficient 
of  order  1019  cm2  s-1.  This  is  considerably  smaller  than  the 
estimate  by  Bell  (2004).  Taking  a  characteristic  CR  energy  of 
10 15  eV,  Bell  (2004)  finds  a  typical  growth  time  for  magnetic 
field  of  order  100  yr.  This  would  yield  a  characteristic  length  scale 
for  magnetic  field  amplification  of  ~  10 18  cm  for  a  3000  km  s-1 
shock,  requiring  a  CR  diffusion  coefficient  of  ~3  x  1026  cm2s  \ 
We  suspect  that  our  simple  estimate  reflects  the  growth  rate  while 
the  shock  may  be  considered  quasi-parallel,  and  that  magnetic 
field  amplification  slows  down  considerably  once  it  becomes 
quasi-perpendicular.  Therefore,  in  taking  a  characteristic  CR  en¬ 
ergy  of  10 15  eV,  Bell  (2004)  is  taking  the  lowest  energy  CRs  for 
which  the  shock  may  be  considered  quasi-parallel,  and  this  result 
may  be  considered  more  realistic. 

Further,  in  Ghavamian  et  al.  (2007)  we  argued  that  the  depth 
of  the  CR  precursor  over  which  electron  heating  occurs  could  not 
be  larger  than  ~10 %vslne  cm,  otherwise  neutral  hydrogen  would 
not  survive  to  encounter  the  shock  front.  We  suggest  here  that 
lower  hybrid  waves  accelerate  the  small  fraction  of  electrons  that 
happen  to  be  in  resonance,  and  that  these  accelerated  electrons 
communicate  their  energy  to  the  rest  of  the  thermal  population 


This  is  the  maximum  electric  field  before  ion  trapping  and  heat¬ 
ing  occurs.  Laming  &  Lepri  (2007  and  references  cited  therein) 
demonstrate  that  when  ijjl\f2k\\yte  Lu/v/2k±vti,  ions  are  heated 
more  effectively  than  electrons  above  this  threshold.  In  our  case, 
the  ions  that  are  heated  will  be  the  lower  energy  part  of  the 
suprathermal  ion  distribution  reflected  from  the  shock,  i.e.,  those 
below  the  injection  threshold  for  diffusive  shock  acceleration  in 
equation  (1)  or  any  of  its  modifications  subsequent  to  the  treat¬ 
ment  of  the  reactive  lower  hybrid  wave  instability  in  §  2.2.  With 
the  wave  electric  field  given  by  equation  (27),  the  electron  mo¬ 
mentum  diffusion  coefficient  in  lower  hybrid  turbulence  varies  as 
v2,  yielding  a  constant  degree  of  heating  with  shock  velocity  if 
the  time  spent  in  the  turbulence  varies  as  1/z;2,  which  would  be 
the  case  if  the  CRs  are  obeying  a  diffusion  law. 

4.  SUMMARY 

We  have  considered  in  more  detail  the  speculation  of  Ghavamian 
et  al.  (2007)  that  lower  hybrid  waves  generated  in  a  cosmic-ray 
precursor  could  be  responsible  for  the  electron  heating  at  colli¬ 
sionless  shocks  in  supernova  remnants.  We  find  that  there  do 
exist  growing  modes  for  the  resonant  or  kinetic  case,  and  that 
the  growth  rate  in  this  case  may  be  sufficient  both  to  survive  the 
damping  by  electrons  and  to  compete  with  magnetic  field  am¬ 
plification  by  modified  Alfven  waves.  Below  a  certain  Alfven 
Mach  number  (roughly  estimated  to  be  ~  12-60)  the  lower  hybrid 
wave  growth  rate  exceeds  that  of  the  modified  Alfven  waves.  The 
modified  Alfven  wave  generation  exists  for  all  magnetic  field 
orientations  with  respect  to  the  shock,  but  is  most  effective  for 
quasi-parallel  case  and  always  generates  new  perpendicular 
field.  Lower  hybrid  waves,  on  the  other  hand,  require  quasi¬ 
perpendicular  field  geometry  in  order  to  grow.  Thus,  a  schematic 
picture  emerges  in  which  far  ahead  of  the  high-Mach  number 
shock,  modified  Alfven  waves  generate  perpendicular  field,  re¬ 
ducing  the  effective  Mach  number  closer  to  the  shock  front  and 
thus  allowing  lower  hybrid  wave  growth  to  occur  in  a  short  re¬ 
gion  before  the  shock  and  to  heat  the  resonant  electrons.  A  critical 
Alfven  Mach  number  around  15  suggests  magnetic  field  ampli¬ 
fication  by  about  an  order  of  magnitude,  similar  to  what  a  com¬ 
parison  of  the  surveys  of  Ghavamian  et  al.  (2007)  and  Schwartz 
et  al.  (1988)  would  suggest,  taking  in  both  cases  the  shock  ve¬ 
locity  where  (Te/Tp)0  ~  1  starts  to  break  down  as  that  where  the 
laminar  shock  approximation  ceases  to  hold. 


No.  1,  2008 


THERMAL  ELECTRON  HEATING  IN  COLLISIONLESS  SHOCKS 


355 


We  have  concentrated  on  the  generation  of  lower  hybrid  waves, 
since  for  these  the  group  velocity  can  be  equal  to  the  shock  ve¬ 
locity  itself,  meaning  that  the  waves  can  stay  in  contact  with  the 
shock  for  long  time  intervals  and  in  principle  grow  to  large  am¬ 
plitudes.  However,  other  wave  modes  that  heat  electrons  are 
certainly  possible,  and  these,  such  as  the  Landau  damping  of 
kinetic  Alfven  waves  (e.g.,  Vinas  et  al.  2000),  do  not  require 
perpendicular  shocks  as  lower  hybrid  waves  do.  In  fact,  Bykov 
&  Uvarov  ( 1 999)  studied  the  generic  case  of  heating  by  turbulent 
modes  in  the  shock  precursor  and  did  identify  an  area  of  pa¬ 
rameter  space  for  which  a  near-inverse-square  relationship  be¬ 
tween  ( Te/Tp)0  and  shock  velocity  could  be  accommodated.  Our 
model  requires  that  cosmic-ray  ions  be  essentially  ubiquitous  at 


SNR  shocks,  with  number  densities  estimated  by  various  means 
in  §  3.  In  a  wider  context,  the  idea  that  cosmic  rays  are  respon¬ 
sible  for  electron  heating  at  fast  shocks  reinforces  the  idea  that 
cosmic  rays  are  an  intrinsic  component  of  the  collisionless  shock 
phenomenon. 
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APPENDIX  A 

COSMIC-RAY  DIFFUSION  COEFFICIENTS 

The  parallel  spatial  cosmic-ray  (CR)  diffusion  coefficient  is  most  easily  obtained  from  its  relation  to  the  pitch-angle  scattering  dif¬ 
fusion  coefficient  in  momentum  space.  The  diffusion  coefficient  in  momentum  space  is  expressed  most  generally  as  (Melrose  1986) 


D 


=  £ 


87i-<72  RM(k) 
h  u>M(k) 


•  v(k,  p,s)\26(ujm  —  sQ  —  £p|j)  A2  ApN^k) 


dik 

(2^F 


(Al) 


where  N ^  is  the  number  density  of  wave  quanta,  RM  is  the  ratio  of  electric  energy  to  total  energy  in  the  wave,  such  that 
Rm  J N^huiM  d2k/(2ir)3=  SE2/&n,  e  is  the  wave  polarization  vector,  and  v  is  the  CR  velocity.  For  pitch-angle  scattering  by  parallel- 
propagating  Alfven  waves,  1  =  a,  so 
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With  u>M  =  &||ha,  Rm  =  (v|)/(2c2),  and  e  •  v  =  v±/2, 
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where  we  have  put  s  =  1  and  taken  ^  O. 

We  now  express  D»  in  terms  of  Daa  by  writing 

f(p,a)  =  fo( p)+ fi(p)  cos  a  +  (l/2)/2(p)cos2a+  .  .  . 

and  substituting  into  the  diffusion  equation 
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Upon  integrating  the  result  over  cos  a  we  obtain,  with  vz  =  v  cos  a, 

dfo  1  df2  v  dfi  = 

dt  3  dt  3  dz 

Multiplying  each  side  by  cos  a  and  then  integrating  over  cos  a  yields 
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With /o  »/i  S>/2,/2  —  —  {2v/3)(dfo/dz)  I  cos  a  s'm2aDaa  d(cos  a)  I  ,  which  when  substituted  into  equation  (A6)  allows  the 
identification 
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With  equation  (A3), 
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This  is  a  factor  of  2-7T  larger  than  the  equivalent  expression  given  by  Blandford  &  Eichler  (1987),  due  to  a  different  definition  of  Um- 
Where  Um  <x  k^,D\\  ex  p2~^,  which  evaluates  to  D\\  ex  vpl/3  orD||  ex  vpm  for  Kolmogorov  or  Kraichnan  turbulence,  respectively.  If 
v  ~  c,  the  dependence  of  Z)  y  on  p  can  usually  be  neglected. 

The  perpendicular  spatial  CR  diffusion  coefficient  has  been  given  in  terns  of  D  y  by  various  authors.  Based  on  numerical  experiments, 
Marcowith  et  al.  (2006)  give  D±  =  r]2+tD y,  where  q  =  SB2 /(SB2  +  ( B )2)  and  the  CR  distribution  function f(p)  cx  p~4~e.  Shalchi  & 
Kourakis  (2007)  and  Zank  et  al.  (2006)  give  D j_  ex  (SB2  /B jj)  2/3Z)j'/3  (Imv)2®  from  analytic  considerations,  where  U d  is  the  two-dimensional 
bend-over  length  scale,  the  inverse  of  the  wavenumber  where  the  inertial  range  onsets,  and  consequently  has  an  even  smaller  dependence  on 
the  CR  momentum  than  the  parallel  diffusion  coefficient  for  relativistic  CRs,  and  has  the  same  dependence  in  the  nonrelativistic  case. 

APPENDIX  B 

GROWTH  RATE  FOR  AN  ELECTROMAGNETIC  INSTABILITY 


For  completeness,  we  give  here  a  treatment  of  the  growth  rate  due  to  CRs  of  electromagnetic  waves  with  frequency  in  the  lower 
hybrid  range  and  show  that  it  is  significantly  smaller  than  either  the  electrostatic  instability  or  the  growth  of  modified  Alfven  waves.  It 
is  relatively  easy  to  show  that  the  reactive  instability  of  Bell  (2004)  has  higher  thresholds  and  lower  growth  rates  as  the  frequency  of 
the  electromagnetic  wave  increases  first  above  the  proton  gyrofrequency  and  then  above  the  electron  gyrofrequency.  Here,  we  con¬ 
centrate  on  the  kinetic  instability  that  might  generate  electromagnetic  waves  in  the  lower  hybrid  range,  whistlers,  adapting  the  ex¬ 
pressions  in  Bell  (2004)  and  Achterberg  (1983), 


Evaluating  the  two  CR  terms, 
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explicitly  assuming  a  K-distribution  for  CRs  in  equation  (B3).  After  some  algebraic  manipulation,  an  integration  by  parts,  the  rewriting 
of/as  a  K-distribution,  and  making  the  substitutions  p\  =  P  and  dP  =  2p±dp± ,  the  right-hand  side  of  equation  ( B3)  can  be  written  as 
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The  integral  over  dP  can  be  evaluated  using  a  hypergeometric  function  (Gradshteyn  &  Ryzhik  1965,  eq.  [3.197. 1])  to  give 
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where  5(1 , 1  +  k)  is  the  beta  function.  Considering  only  the  kinetic  case,  using  the  Landau  prescription  for  this  integral  with  the  pole 
at  u>  —  k ||ti||,  the  imaginary  part  (i.e.,  the  portion  relevant  for  the  growth  rate)  is 
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Setting  ui  — >  u>  +  rys  in  the  dispersion  relation  and  taking  only  the  imaginary  parts,  we  get 


no  ,  o  ,  Mat  3  2,  (~ium 

0  =  lli'Yg  +  2urfg  ±  —7s 7777  +  (2tt)  <7  k\\  I  — - mv. 


to2  'y  neq  '  1  11 V  ^|| 


(2kP2) 


(2  Kp2  +  m~v 


?  2 

(2k  -  3)T(k) 


^tB(  1,1+k) 


v\  _)_m_  (  _  2kP2  +  m2v2\  nc  r _ 

c2  m2  k\\2  \  ’  ™2c2  )  4x/2(7tk)3/2;?3  T(k  -  1/2) 


(B7) 


assuming  7  w/A:||  <SC  and  including  the  normalization  of /  and  the  factor  of  k  y  /2  that  were  omitted  during  the  evaluation  of  the  integral. 

From  this  we  have 


7s-  “ 


'/2ncR  fl?  ( yuim/ka  —  mvs ) 

«,  2w  Pt 


2k  -  3  T(k) 

k>/2(k+  1)T(k-  1/2) 


2F1(2,1;2  +  k;1). 


(B8) 


Electromagnetic  waves  in  the  lower  hybrid  frequency  range  are  parallel-propagating  whistlers,  with 


2A 


'2  uj{kle  —  ui)  ’ 
2  uj/k 


1  +  k2c2/u)2e  ' 

We  assume  dui/dk  oc  cL>/dA:||  ~  us,  and  hence,  w/A:||  ~  vs/2  for  A:  <S  c!u>pe.  Thus,  for  k  =  2 

2«cr  it2 


7s-' 


3k,-  w 


(>-!)■ 


(B9) 

(BIO) 


(Bll) 


This  is  significantly  smaller  than  the  growth  of  lower  hybrid  waves,  which  is  of  order  H,;7cr/k;-.  Further,  since  whistlers  carry  energy 
along  magnetic  field  lines,  like  Alfven  waves,  only  for  specific  shock  obliquities  will  the  energy  of  the  waves  stay  in  contact  with  the 
shock  and  allow  large  wave  intensities  to  build  up.  Electromagnetic  waves  with  frequency  above  the  electron  gyrofrequency  ( O -  and 
X-modes)  have  phase  velocities  greater  than  c,  and  so  cannot  be  excited  by  kinetic  instabilities. 
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